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Abstract: We study dynamical effects of introducing noncommutativity on string 
worldsheets by using a matrix model obtained from the zero-volume limit of 
four-dimensional SU(iV) Yang-Mills theory. Although the dimensionless noncom- 
^ ■ mutativity parameter is of order 1/N, its effect is found to be non-negligible even 

in the large N limit due to the existence of higher Fourier modes. We find that the 
Poisson bracket grows much faster than the Moyal bracket as we increase N, which 
means in particular that the two quantities do not coincide in the large N limit. The 
well-known instability of bosonic worldsheets due to long spikes is shown to be cured 
by the noncommutativity. The extrinsic geometry of the worldsheet is described by 
a crumpled surface with a large Hausdorff dimension. 
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1. Introduction 

It is well known that the traditional bosonic worldsheet theories described for exam- 
ple by the Nambu-Goto action, the Polyakov action and the Schild action are not 
well-defined non-perturbatively, as was first observed in the large dimension limit by 
Alvarez [1]. This conclusion has been confirmed more directly by the dynamical tri- 
angulation approach [2, 3], where the discretized worldsheet embedded in the target 
space degenerates to long spikes, and hence one cannot view it as a proper approxi- 
mation of a continuous worldsheet. It is therefore of interest to inquire whether it is 
possible to have different types of string models, where the worldsheet is well-defined 
in the sense that it does not degenerate into long spikes. Some attempts have been 
made to make the string more "stiff" by the introduction of extrinsic curvature [4], 
but it is not clear whether this will work non-perturbatively [5]. 

In the present paper we have investigated whether the introduction of noncom- 
mutativity on the worldsheet changes the situation. We study such a system by using 
a matrix model with N x N hermitian matrices with pi — 1, 2, 3, 4. Corresponding 
to each matrix one can construct a field A M (<r) by the Weyl transformation. Here 
a stands for the two discrete variables o\ and 02, representing the worldsheet coor- 
dinates. Then there exists the following relation between the partition functions in 
the matrix model and in the corresponding Weyl-transformed model, 

J d^exp (_L T r[^, A,] 2 ) = J VX,(a)e~ s , (1.1) 
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where the action S for the latter is given by 

S = -^nY,{ X M* X M - X u {a)*X,{<y)) 2 . (1.2) 

Here * denotes the star product with the dimensionless noncommutative parameter 
being of order 1/JV. The partition function on the left hand side of eq. (1.1) was 
shown to be finite for N = 3,4, 5 numerically by Ref. [6]. Monte Carlo simulations 
up to N = 256 [7] further suggests that the statement extends 1 to iV = oo. This 
means that the corresponding Weyl-transformed model is also well-defined in the 
large N limit. 

The partition function on the right hand side of eq. (1.1) defines a noncommu- 
tative two-dimensional field theory, which is invariant under the star-unitary trans- 
formation, 

X tl {a)^g*{a)^X,{a)^g{a) , with g*(a) * g(a) = g(a) * g*(a) = 1 . (1.3) 
The action (1.2) resembles the Schild action 

S s hid 1 /dV ( dX » dX » _ dx » dx t> 
g 2 N J \ dai da 2 do\ da 2 

the difference being that the Poisson bracket has been replaced by the Moyal bracket 
in the noncommutative case. It has been pointed out long time ago [9] that if the 
higher modes in the mode expansion of X^(a) can be neglected, the star-Schild 
action in (1.2) reduces to the Schild action (1.4) in the large N limit. The issue is 
also relevant when one considers matrix models as a nonperturbative definition of 
M-theory and type IIB superstring theory [10, 11]. Of course, whether the higher 
modes can really be neglected or not is a non-trivial dynamical question. One of the 
main results of this paper is that for the bosonic model we consider, the Poisson and 
the Moyal brackets are very different due to the effect of the higher modes, and that 
this effect increases with increasing N. To our knowledge, this is the first time that 
the agreement (or disagreement) of the two actions in the large N limit is discussed 
in a dynamical context. 

Recently the Schild action has been investigated non-perturbatively from the 
point of view of dynamical triangulation [3]. The result is similar to the one ob- 
tained for the Nambu-Goto action, namely that it is dynamically favorable to have 
a "surface" which degenerates to long spikes, and hence the notion of a worldsheet 
looses its meaning 2 . 

1 This has been proved rigorously in a preprint [8], which appeared while our paper was being 
revised. 

2 A similar result is not valid in the supersymmetric case, where the worldsheet exists if the 
fermions are coupled to a bosonic Schild action [3]. However, in the present paper we shall discuss 
the bosonic case only. 
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We shall take a different approach and ask whether the star-Schild action in 
(1.2) provides a bosonic string theory with a well-defined worldsheet. The effects of 
noncommutativity are found to be drastic. The average link length is finite and it is 
observed to be considerably smaller than the average extent of the surface, which is in 
sharp contrast [1, 3] to the conclusion for the Schild action (1.4) that the worldsheet 
becomes completely unstable due to long spikes. This is not a contradiction since 
the Poisson and Moyal brackets are indeed found to be quite different for the theory 
under consideration. On the other hand, the extrinsic geometry of the worldsheet 
for the star-Schild action is described by a crumpled surface with a large Hausdorff 
dimension. 

The matrix model on the left hand side of eq. (1.1) can be regarded as the 
hermitian matrix version of the Eguchi-Kawai model [12]. (See Refs. [13] and [14] 
for the hermitian matrix version of the Eguchi-Kawai model with quenching [15] or 
twists [16], respectively.) Large N scaling of SU(iV) invariant correlation functions 
was obtained analytically to all orders of the l/D expansion [7], and the results were 
reproduced by Monte Carlo simulations [7, 17]. In particular, the one-point function 
of a Wilson loop was observed to obey the area law, which suggested that the model 
is actually equivalent to large N gauge theory for a finite region of scale even without 
quenching [15] or twists [16]. The noncommutative worldsheet theory studied in the 
present paper is therefore related to the large N QCD string (For a recent review, 
see [18].). 

The organization of this paper is as follows. In Section 2 we briefly review the 
map from matrices to noncommutative fields. In Section 3 we interpret the matrix 
model as noncommutative worldsheet theory. We then discuss the star-unitary in- 
variance and the important question of gauge fixing. In Section 4 the results are 
presented. Section 5 contains the summary and discussions. In an Appendix we 
make some comments on the relationship between the star-unitary invariance and 
the area-preserving diffeomorphisms. 

2. From matrices to noncommutative fields 

In order to derive the equivalence (1.1) between the matrix model and the noncom- 
mutative worldsheet theory, we briefly review the one-to-one correspondence between 
matrices and noncommutative fields. Most of the results in this Section are known 
in the literature (see e.g., [19, 20]), but are given in order to fix the notation. 

Throughout this paper, we assume that N is odd 3 . We first introduce the 't 

3 We expect that the large TV dynamics of the matrix model on the l.h.s. of (1.1) is independent 
of whether N is even or odd. This has been also checked numerically for various SU(7V) invariant 
quantities. However, the one-to-one correspondence between the matrices and noncommutative 
fields works rigorously only for odd N. 
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Hooft-Weyl algebra 

r 1 r 2 = cur 2 r 1 , (2.1) 

where to = e Am / N . It is known that the representation of Tj using N x N unitary 
matrices is unique up to unitary transformation Tj — > UFiW, where U G SU(iV). 
An explicit representation can be given by 



/0 1 

1 



Vi 



1 

0/ 



r. 



V 



ur 



\ 



(N-l) J 



(2.2) 



(2.3) 



Then we construct N x N unitary matrices 

J k = (r 1 ) ki (r 2 ) ka e ~^ ik ^/ N , 

where A; is a 2d integer vector and the phase factor e ~ 2mk i k 2/ N [ s included so that 

J-k = {JkY ■ (2.4) 
Since (Ti) N = 1, the matrix Jk is periodic with respect to fcj with period N, 

Jki+N,k 2 = Jki,k 2 +N — Jkifa ■ (2-5) 

We then define the N x N matrices 

A(<t) = Jk ' ( 2 - 6 ) 

ki&Z N 

labelled by a 2d integer vector a. Note that A(cr) is hermitian due to the property 
(2.4). It is periodic with respect to ui with period N, 



A(oi + N, a 2 ) = A(oi, a 2 + N) = A(a u a 2 ) . 
It is easy to check that A(er) possesses the following properties. 



tr Ai 



N 

N 2 1 N 



1 



tr(A(<x)A(<7')) =7V 2 ^ (mod 



JV) • 



(2.7) 



(2.8) 
(2.9) 

(2.10) 



Now let us consider gl(iV,C), the linear space oi N x N complex matrices. The 
inner product can be defined by tr(M|M 2 ), where M 1 ,M 2 e gl(iV,C). Eq. (2.10) 
implies that A(cr) are mutually orthogonal and hence linearly independent. Since 
the dimension of the linear space gl(iV,C) is iV 2 and there are iV 2 A(a)'s that are 
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linearly independent, A(cr) actually forms an orthogonal basis of gl(N, C). Namely, 
any N x N complex matrix M can be written as 

m = jp E /(^ A (°-) ? (2-11) 

where /(er) is a complex- valued function on the 2d discretized torus obeying periodic 
boundary conditions. Using orthogonality (2.10), one can invert (2.11) as 

/(a) = itr(MA(,)) . (2.12) 

The fact that (2.11) and (2.12) hold for an arbitrary N x N complex matrix M 
implies that A(er) satisfies 

± Ay((7)AH((7) = NS u S jk . (2.13) 



Note also that using (2.8) with (2.11) or using (2.9) with (2.12), one obtains 

— tr M = — 

N N 2 



■^i ^ E /(*)■ ( 2 - 14 ) 



We define the "star-product" of two functions /i(<r) and /2(c) by 

fi{°) * /2W = ^ tr(M 1 M 2 A(cr)) , (2.15) 

where 

M « = E /<» A ( (J ) « = 1, 2 . (2.16) 
The star-product can be written explicitly in terms of f a (&) as 

hi?) * / 2 W = ^ E E /i(0/2(^)e- 2 ^ fc ^-^ )fe - CTfc)/JV , (2.17) 

where is an antisymmetric tensor with 612 = 1. This formula can be derived in 
the following way. Substituting (2.16) into (2.15) and using the definition (2.6) of 
A(cr), one obtains 

From the definition of J k , one easily finds that 

1 tr(J p J q J k ) = e 2m ^^ N 5 k+p+qfi(mod N) . (2.19) 
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Integration over k, p and q in (2.18) yields eq. (2.17). 

In order to confirm that the star-product (2.17) is a proper discretized version 
of the usual star-product in the continuum, we rewrite it in terms of Fourier modes. 
We make a Fourier mode expansion of f a (cr) as 

/«(*) = E /«(*) e2mk ' a/N ( 2 - 2 °) 
/«(*) = ^ E /» ^ 2mk ' a,N ■ ( 2 - 21 ) 

Integrating (2.18) over k, £ and 77, one obtains 

fi(<r)*f 2 (a)= E E /i(P)/a(9) e**^/^^^ , (2.22) 

which can be compared to the usual definition of the star-product in the continuum 
/i(<t)*/ 2 (<7) = fiWexpfyei&dfjMa) . (2.23) 

In the present case, the noncommutativity parameter is of order 1/N, and therefore 
the star-product reduces to the ordinary product in the large N limit if / a (c) contains 
only lower Fourier modes (p^, qj <C \/~N). 

It is obvious from the definition (2.15) that the algebraic properties of the star- 
product are exactly those of the matrix product. Namely, it is associative but not 
commutative. Note also that due to (2.14), summing a function f(a) over a corre- 
sponds to taking the trace of the corresponding matrix M. Therefore, 

E /i(*)*/2(<r)*---*/n(<7) (2.24) 

is invariant under cyclic permutations of / Q (cr). What is not obvious solely from the 
algebraic properties is that 

E /i(') */:»(*) = E » ( 2 - 25 ) 

which can be shown by using the definition (2.15) with eq. (2.9). 
For later convenience, let us define the Moyal bracket by 

/ 2 (*)1 = i ^ * /2W - /2(a) * A(a)) 

= E /i(P)/»(9) sin (^f^) e 2 -^)'^ . (2.26) 

P9 
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We also define the Poisson bracket on a discretized worldsheet. Namely when we 
define the Poisson bracket 



UMJii*)} = e i3 , (2.27) 



'3 

we assume that the derivatives are given by the lattice derivatives 

^ = +«)-/(-!)). (2-28) 

where a = 2n/N is the lattice spacing. The Poisson bracket thus defined can be 
written in terms of Fourier modes as 

{fMJM} = - (^) 2 £ /i(p)/2(?h-*sin (^) sin (^) e-^W- . 



(2.29) 

Note that the appearance of sines is due to discretization of the worldsheet. The 
Moyal bracket (2.26) and the Poisson bracket (2.29) agree in the large N limit if 
nonvanishing Fourier modes are those with pj, qj <C \fN. 



3. Matrix model as noncommutative worldsheet theory 

Let us proceed to the derivation of the equivalence (1.1) between the matrix model 
and the noncommutative worldsheet theory. We start from the matrix model with 
the action 

S = --^tr([A lt ,A v ] 2 ), (3.1) 

which can be obtained from the zero-volume limit of D-dimensional SU(iV) Yang- 
Mills theory 4 . The indices fi run from 1 to D. 

As we have done in (2.11), we write the N x N hermitian matrices in (3.1) 

as 

A " = Jp E X » A (^) > ( 3 - 2 ) 

where X^{a) is a field on the discretized 2d torus obeying periodic boundary condi- 
tions. Since is hermitian, X^a) is real, due to the hermiticity of A (a). Eq. (3.2) 
can be inverted as 

X,(a) = -tr(A,A(a)) . (3.3) 

We regard a as (discretized) worldsheet coordinates and X^{a) as the embedding 
function of the worldsheet into the target space. 

4 For D = 10, the action (3.1) is just the bosonic part of the IIB matrix model [11]. The 
dynamical aspects of this kind of matrix models for various D with or without supersymmetry have 
been studied by many authors [6, 7, 17, 21] both numerically and analytically. 
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Using the map discussed in the previous section we can rewrite the action (3.1) 
in terms of X^cr) as 5 

= ^(|V (3-4) 

If X m (<t) are sufficiently smooth functions of a, the Moyal bracket can be replaced 
by the Poisson bracket, and therefore we obtain the Schild action 

Sschiid = i(§V E {X^lXMY ■ (3-5) 

Let us discuss the symmetry of the theory (3.4), which shall be important in 
our analysis. For that we recall that the model (3.1) is invariant under the SU(iV) 
transformation 

A^gA^gK (3.6) 

From the matrix-field correspondence described in the previous section, one easily 
finds that the field X^a) defined through (3.3) transforms as 

X^a)^g(a)^X^a)*g(aY , (3.7) 

where g(a) is defined by 

g(tr) = ± ti(g . (3.8) 

The fact that g e SU(iV) implies that g(a) is star-unitary; 

g{a)**g(a)=g(a)*g{ay = l. (3.9) 

The action (3.4) is invariant under the star-unitary transformation (3.7) as it should. 
We shall refer to this invariance as 'gauge' degrees of freedom in what follows. 

Even if X^a) is a smooth function of a for a particular choice of gauge, it can 
be made rough by making a rough star-unitary transformation. Let us quote an 
analogous situation in lattice gauge theory. In the weak coupling limit, the configu- 
rations can be made very smooth by a proper choice of the gauge. However, without 
gauge fixing, they are as rough as could be due to the unconstrained gauge degrees 
of freedom. Similarly when we discuss the smoothness of X M (cr), we should subtract 
the roughness due to the gauge degrees of freedom appropriately. Therefore, a natu- 
ral question one should ask is whether there exists at all a gauge choice that makes 
Xn(a) relatively smooth functions of a. 

5 While this work was being completed, we received a preprint [22] where the worldsheet theory 
(3.4) is discussed from a different point of view. 
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Figure 1: A typical TV = 35 X^(a) (// = 1) configuration after the gauge fixing. 



In order to address this question, we specify a gauge-fixing condition by first 
defining the roughness of the worldsheet configuration X^a) and then choosing a 
gauge so that the roughness is minimized. A natural definition of roughness is 



1 

2PP 



o-ieZiv i =1 



(3.10) 



which is Lorentz invariant. The gauge fixing is analogous to the Landau gauge in 
gauge theories. The roughness functional (3.10) can be written conveniently in terms 
of as (See Appendix A for the derivation.) 



-E 

2N ^ 

u 



4 sin' 



tt(I-J) 



N 



\{A^) U \ 2 + (A^u - (A^) I+ N-L J+ N^l 



(3.11) 



4. Results 



Our numerical calculation starts with generating configurations of the model (3.1) 
for D = 4 and N = 15, 25, 35 using the method described in Ref. [7]. For each config- 
uration we minimize the roughness functional / defined by eq. (3.11) with respect to 
the SU(A^) transformation (3.6). We perform 2000 sweeps per configuration, where a 
sweep is the minimization of / with respect to all SU(2) subgroups of SU(A r ) [7, 17]. 
From the configuration A^ obtained after the SU(A^) transformation that minimizes 
/, we calculate through (3.3) the worldsheet configuration X^(a). When we define 
an ensemble average ( • ) in what follows, we assume that it is taken with respect 
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Figure 2: The effect of cutting off Fourier modes higher than k c on the worldsheet of 
Fig. 1. On the left k c = 1, on the right k c = 3 and on the bottom k c = 5. 

to Xfj,(a) after the gauge fixing. The number of configurations used for an ensemble 
average is 658, 100 and 320 for N = 15, 25, 35 respectively. We note that in the 
present model, finite N effects is known [7] to appear as a 1/N 2 expansion. Also, the 
large N factorization is clearly observed for N = 16,32 [7]. (We have checked that 
it occurs for iV = 15, 25, 35 as well.) We therefore consider that the N we use in the 
present work is sufficiently large to discuss the large N asymptotics. 

We also note that the parameter g, which appears in the action (3.4), can be 
absorbed by the field redefinition X'^a) = ^=X^(cr). Therefore, g is merely a scale 
parameter, and one can determine the g dependence of all the observables on dimen- 
sional grounds. The results will be stated in such a way that they do not depend on 
the choice of g. 

In Fig. 1 we show a typical worldsheet configuration for N = 35 after the gauge 
fixing. We observe that the worldsheet has no spikes. We compute the Fourier modes 
X^(k) of X^(a) through 

= E M^- 2 ^ = ^ tr(A„ J k ) , (4.1) 

where the range of k\ and A; 2 are chosen to be — (JV — l)/2 to (N — l)/2. Fig. 2 
describes how the same worldsheet configuration shown in Fig. 1 looks like if we cut 
off 6 the Fourier modes higher than k c . We find that the configuration obtained by 

6 More precisely, we keep the modes X^k) with k\ < k c and k% < k c and set the other modes to 
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15 25 35 

N 

Figure 3: The fluctuation of the surface R 2 and the roughness functional /, which repre- 
sents the average link length in the target space, are plotted against N with the normaliza- 
tion factor g~ x . The solid straight line is a fit to the power- law behavior (R 2 ) oc 5 ArO-493(3) _ 
The dotted line is drawn to guide the eye. 



keeping only a few lower Fourier modes already captures the characteristic behavior 
of the original configuration. We have checked that this is not the case if we do not 
fix the gauge. 

We measure the fluctuation of the surface, which is given by 7 

The average of the fluctuation is finite for finite N, and it is plotted in Fig. 3. The 
power-law fit to the large N behavior (R 2 ) ~ gN p yields p = 0.493(3), which is 
consistent with the result p — 1/2 obtained in Ref. [7]. Although the fmiteness of 
the fluctuation already implies a certain stability of the worldsheet, we note that the 
fluctuation defined by the l.h.s. of (4.2) is actually invariant under the star-unitary 
transformation (3.7). In particular, the fluctuation is finite even before the gauge 

zero by hand. 

7 We note that the action (3.1) is invariant under i— > + a^ljv- Therefore, the trace part 
of A^ is completely decoupled from the dynamics. We fix these degrees of freedom by imposing 
the matrices A^ to be traceless. In the language of the worldsheet theory (3.4), the symmetry 
corresponds to the translational invariance X fl (a) X^(a) + const.. The tracelessness condition 
for A,, maps to £ CT . eZjv X^(a) = 0. 
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Figure 4: The Fourier mode amplitudes (^X^(k)X^(—k)^ /(gN 1 / 2 ) are plotted against k\ 
(k 2 = ki) for TV = 15,25,35. 



fixing. Therefore, the smoothness of X^ia) (which appears only after the gauge 
fixing) is a notion which is stronger than the finiteness of the fluctuation. 

Let us point out also that the roughness functional I actually represents the 
average link length in the target space. We therefore plot (I) in Fig. 3 and compare 
it with the fluctuation (4.2). The former is observed to be smaller than the latter 8 
which is consistent with the observed smoothness of X^a). The ratio (I) / (R 2 ) is 
0.364(1), 0.338(2), 0.3290(5) for N = 15,25,35, respectively, which may be fitted 
to a power-law behavior (/) / (R 2 ) ~ jSf-°- 120 ( 4 ) . This indicates a tendency that the 
worldsheet is getting smoother as N increases. 

In order to quantify the smoothness of X^a) further, let us examine the Fourier 
mode amplitudes. We first note that there is a relation 

(x,(k)X,(-k)) = I-L J2 X»A ~ 0(^/2) . (4 . 3) 

fcieZjv \ o-iGZjv / 

This motivates us to plot <^X M (A;)X M (— k)^ /(gN 1 / 2 ) with this particular normaliza- 
tion. The results are shown in Fig. 4 and Fig. 5. We find a good scaling in N; 
data points for different iV fall on top of each other. The discrepancy in the large k 
region can be understood as a finite iV effect. The k dependence of the Fourier mode 

we do not fix the gauge, we observe that the two quantities (R 2 ) and (7) coincide, meaning 
that the X fl (<r) is completely rough. 
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Figure 5: The Fourier mode amplitudes (x^tyX^—k)^ / (gN 1 / 2 ) in Fig. 4 are now- 
plotted in the log-log scale in order to visualize the power-law behavior. The straight line 
is a fit to C\k\~i with q = 1.96(5) for the N = 35 data. 



amplitudes suggests that the higher modes are indeed suppressed. Moreover we find 
that there exist a power-law behavior 9 

(x^X.i-k)) ~ const. ■ \k\~« , (4.4) 

where \k\ — a/ (&i) 2 + (^) 2 - Assuming that the constant coefficient on the r.h.s. of 
(4.4) is independent of N, as suggested by the observed scaling in N, the power q 
must be q > 2 in order that the sum on the l.h.s. of (4.3) may be convergent in 
the large N limit. The power q extracted from N = 35 data is q = 1.96(5), which 
may imply that we have not reached sufficiently large k (due to the finite iV effect 
mentioned above) to extract the correct power. Although we have seen that the 
amplitudes of the higher Fourier modes are suppressed, we should remember that 
their number for fixed \k\ grows linearly with \k\. Therefore, the higher modes can 
still be non-negligible. 

Let us turn to the question whether the action S in (3.4) approaches the Schild 
action S's c hiid m (3-5) in the large iV limit. In terms of Fourier modes, the two actions 

9 If we do not fix the gauge, we observe that the r.h.s. of (4.4) is replaced by a constant indepen- 
dent of fc, which means that the worldsheet behaves like a white noise. The constant is proportional 
to 1/N 2 , as expected from the relation (4.3), which is gauge independent. Note, in particular, that 
the scaling behavior observed in Fig. 5 emerges only after the gauge fixing. 



13 



100000 



10000 



1000 



100 





<s> 




<S Schild > °" 


7 

















15 



25 



35 



N 



Figure 6: Plot of (S) and (S'schild) as a function of N. The solid line represents the exact 
result (S) = (N 2 — 1). The dotted line is drawn to guide the eye. 
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^rl^j e rs sin ^Pr^ sin (jf( k + 1 + ( 4 - 6 ) 



xea sin ( —A;,; 1 sin 



If the higher Fourier modes can be neglected one can see from eqs. (4.5) and (4.6) that 
S = S'schild i n the large iV limit. We measure both quantities and plot the results in 
Fig. 6. The average of the action (3.4) is known [7] analytically (S) = N 2 — 1, which 
is clearly reproduced from our data. On the other hand, (Sschiid) is much larger, and 
moreover it grows much faster, the growth being close to 0(iV 4 ). Therefore we can 
safely conclude that the two quantities do not coincide in the large iV limit. 

The disagreement of the two actions (3.4) and (3.5) in the large N limit implies 
that the higher modes play a crucial role. In order to see their effects explicitly, 
we cut off the Fourier modes higher than k c . In Fig. 7, we plot the average of the 
actions thus calculated against k c for N = 35. The two actions with the cutoff at 
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Figure 7: The two actions (Sschiid) ano - (^) are computed by cutting off the Fourier modes 
higher than k > k c by hand. The results are plotted as a function of the mode cutoff k c 
for N = 35. 



k c are indeed identical for small k c , but they start to deviate from each other as k c 
increases, ending up with totally different values at k c = (N — l)/2, i.e. when all the 
modes are included. 

Finally, let us discuss the extrinsic geometry of the worldsheet. One may define 
the Hausdorff dimension d H of the worldsheet through 

^(VM)" as^oo, (4,, 

where R, which is defined by eq. (4.2), represents the extent of the worldsheet in 
the target space. We have defined A, the total area of the worldsheet in the target 
space, by 

A = a 2 J2 \j{X»{v\X u {a)Y , (4.8) 

where the Poisson bracket { • } is defined by (2.27) and a = ^ is the lattice spacing 
on the worldsheet. Eq. (4.8) is nothing but the Nambu-Goto action. The scale 
parameter £ in eq. (4.7) should be introduced in order to make the equation consistent 
dimensionally. A natural choice of the fundamental scale I is the average link length 
in the target space, which is represented by yj (I), where / is the roughness functional 
defined by eq. (3.11). We observe, as expected, that (A) ~ N 2 (I), which means 



15 



that the l.h.s. of (4.7) is of 0(N 2 ). As we have seen in Fig. 3, we observe that 
(I)/(R 2 ) ~ N-°- 120 ^. This means that the Hausdorff dimension du defined by 
(4.7) with the choice i = is d H ~ 33, which might suggest that actually 

d H — oo. Therefore, the extrinsic geometry of the embedded worldsheet is described 
by a crumpled surface. 

5. Conclusions 

In this paper we have studied the star-Schild action (1.2) resulting from the zero- 
volume limit of SU(iV) Yang-Mills theory. We find that the star-Schild action does 
not approach the Schild action (1.4) due to the important role played by the ever- 
increasing number of higher modes. This has some implications to the ideas presented 
long ago [9] that it might be that the Schild action represents the large N QCD 
action. From our results the two actions differ more and more with increasing N . 
The Poisson bracket increases much faster with iV than the corresponding Moyal 
bracket. Our conclusion is therefore that QCD strings would be described by a 
noncommutative string theory defined by the star-Schild action, rather than the 
standard Schild action. 

As we have seen, it is possible to find a star-unitary transformation g{a) such 
that the surfaces defined by the star-Schild action are regular (i.e., do not have long 
spikes), and this action therefore defines a new type of string theory. The theory 
is invariant under star-unitary transformations, which generalize the area-preserving 
diffeomorphisms, the invariance of the usual Schild action. As we discuss in the 
Appendix B, the reparametrization invariance of the worldsheet fields X^a) is re- 
stricted to linear transformations of the cr's in the case of the star-Schild action. 
(Note, however, that the reparametrization invariance of the usual Schild action is 
also much reduced relative to the Nambu-Goto action.) The star-unitary transforma- 
tions transform the worldsheet configuration in such a way that the changes cannot 
be absorbed by a reparametrization. In obtaining a regular surface we have chosen 
a particular "gauge" . The regularity will not be changed drastically under smooth 
star-unitary transformations. However, if we do not fix the gauge, we obtain spiky 
surfaces, which is connected to a regular surface by a rough star-unitary transforma- 
tion g(a). The main point is that it is at all possible to obtain a regular surface by 
fixing the gauge properly. 

A possible intuitive understanding of the regularity of the worldsheet in the 
noncommutative string is that the action contains higher derivatives in the star- 
product. Let us recall that one of the motivations for introducing extrinsic curvature 
(which also contains higher derivatives in a different combination) was [4] that this 
extra term in the action makes the worldsheet more stiff. One would also expect 
a similar effect from the introduction of higher derivatives, because an extremely 
rough worldsheet with long spikes would have at least some derivatives rather large. 
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Although the star-product contains these derivatives in a special combination, it is 
difficult for all these large derivatives to cancel, and hence a surface dominated by 
long spikes would not be preferred. This is confirmed by the observation that the 
average link length is much smaller than the average extent. The extrinsic geometry 
of the embedded surface, on the other hand, is described by a crumpled surface with 
a large Hausdorff dimension. 

It would be of interest to address the issues studied in this paper in the su- 
persymmetric case using the numerical method developed in Ref. [17]. We hope to 
report on it in a future publication. 
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A. Derivation of the roughness functional 

In this appendix, we rewrite the roughness functional (3.10) in terms of the matrices 
and derive (3.11). For that purpose, we introduce N x N unitary matrices 



JV-l 



Di = (rj)T- (A.i) 



JV-l 



which satisfy 

One can check that 

which implies 



D 2 = (T 1 )— , (A.2) 

DjViD) = e- 27vi5 ^ /N ri . (A.3) 

DjA(a)D}^A(a-j) , (A.4) 

X^a + j) = 1 tr(£>^£>jA((7)) . (A.5) 

Thus the matrix Dj plays the role of a shift operator in the j-direction. Now we can 

rewrite the roughness functional I in terms of the matrices as 

3 

= ^££|p;4^5)«-(40«| 2 • (a.6) 

U 3 
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Using the explicit form of T M given by eq. (2.2), we obtain 

(D 1 A,D\) IJ = {AJueW-W* (A.7) 
{D 2 A,D\)jj = (AJ I+EfltJ+Efl , (A.8) 

which yields (3.11). 

B. Star-unitary invariance and area-preserving diffeomorphisms 

In this appendix, we discuss the relationship between the symmetry of the Schild 
action and that of the star-Schild action. Here only we consider that the worldsheet 
is given by an infinite two-dimensional flat space parametrized by the continuous 
variables o\ and a 2 . Let us define the Schild and star-Schild actions 



h = J d 2 a{Mv)M*)} 2 (B-l) 
h = J d 2 v{M<r),M°)}9, (B.2) 



where the Poisson and Moyal brackets are defined by 

Wl(ff) '^ )} = e '^| (R3) 
{{0x(a), M°)}e = ^iW sin(0e^d*)0 2 (<7) . (B.4) 

The Schild action I± is invariant under the area-preserving diffeomorphism 

(Ti^ (n + eijdjf{a) , (B.5) 

where f{a) is some infinitesimal real function of a. Under the infinitesimal area- 
preserving diffeomorphism, the fields transform as a scalar 0q,(<t) = (p' a {a r ), so that 
one can state the invariance as the one under the field transformation 

0» (f> a {a) + {f(a), 0»} . (B.6) 

On the other hand, the star-Schild action I2 is invariant under the star-unitary 
transformation 

0» - <t> a {<j) + If {a), M°)}e ■ (B.7) 

Obviously, this transformation (B.7) reduces to (B.6) if a (a) and f(a) do not contain 
higher Fourier modes compared with Q~ x l 2 . 

In general the two transformations (B.6) and (B.7) differ. However, we note that 
they become identical if f[a) contains terms only up to quadratic in a as 

f(a) = diOi + bijOiOj , (B.8) 
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where is a real symmetric tensor. From (B.5), one finds that the corresponding 
coordinate transformation is 

0i h-> cfi + (vi + \ij<Jj) , (B.9) 

where V; L = ena\ and = eubij, which is traceless. This transformation includes 
the Euclidean group, namely translation and rotation. Thus we find that the linear 
(finite) transformation of the coordinates a\ = AijUj + Vi, where det A = 1, can be 
expressed as a star-unitary transformation. In other words, the reparametrization 
invariance of the star-Schild action is restricted to such linear transformations. 
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